Compressibility and entropy of cold fermions in one dimensional optical lattices 
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We calculate several thermodynamic quantities for repulsively interacting one dimensional 
fermions. We solve the Hubbard model at both zero and finite temperature using the Bethe ansatz 
method. For arbitrary values of the chemical potential, we calculate the particle number density, 
the double occupancy, various compressibilities, and the entropy as a function of temperature and 
interaction. We find that these thermodynamic quantities show a characteristic behavior so that 
measurements of these quantities can be used as a detection of temperature, the metal-insulator 
transition, and metallic and insulating phases in the trap environment. Further, we discuss an 
experimental scheme to extract these thermodynamic quantities from the column density profiles. 
The entropy and the compressibility of the entire trapped atomic cloud also reveal characteristic 
features indicating whether insulating and/or metallic phases coexist in the trap. 



I. INTRODUCTION 

Based on Bethe ansatz solutions of the spin-half 
Heisenberg chain Lieb and Wu proposed solutions 
of the one diniensional(lD) Hubbard model as a set of 
algebraic equations, known as the Lieb-Wu equations Q. 
Since then, the ID Hubbard model has been considered 
as an exactly solvable model and is used as a toy model 
for non-perturbative strongly correlated effects in inter- 
acting electron systems. Later, Takahashi developed a 
classification of the Lieb-Wu equations and introduced a 
set of non-linear integral equations, known as the ther- 
modynamic Bethe ansatz (TBA), to replace the Lieb-Wu 
equations. Any thermodynamic quantity that relates to 
the energy spectrum of the ID Hubbard model can be 
calculated from numerical solutions of the TBA equa- 
tions Ji] . Due to the advent of Feshbach resonances and 
optical lattices in the field of ultra-cold atoms the 
ID Hubbard model has evolved from a toy model to a 
paradigm of experimental relevance for strongly corre- 
lated particle systems. The high degree of tunability and 
the absence of disorder in optical lattice experiments offer 
an opportunity to explore fundamental physical phenom- 
ena in a controlled manner. 

The first successful experimental emulation of the Hub- 
bard model in optical lattices was achieved by Greiner et 
al. [sj using bosonic atoms. This was after the theoretical 
prediction of a Bose-superfiuid to Mott-insulator quan- 
tum phase transition by Jaksch et al Q. With recent 
extraordinary experimental developments in manipulat- 
ing Fermi atoms in optical lattices [rl-fioj. experimental- 
ists are now addressing the rich many-body scenario in 
these systems. Fermions in optical lattices can open new 
regimes and routes for understanding many body physics 
beyond condensed matter models. For example, mixtures 
of two-component atoms with different masses introduce 
spin-dependent hopping amplitudes. This may effect the 
stability of the possible quantum phases or even induce 
new phases. 

Parallel to the experimental progress of manipulating 



fermions in optical lattices, recently there have been the- 
oretical studies on the role of an external harmonic po- 
tential and its effect on compressibility and double occu- 
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The underlying harmonic trapping po- 
tential present in all cold-gas experiments causes density 
to vary across the lattice. This trapping-potential in- 
duced inhomogeneity allows for the spatial coexistence 
of metallic and insulating phases. Measurements of lo- 
cal and/or global thermodynamic properties allows one 
to detect these different phases. In the present study, we 
use an exactly solvable ID Hubbard model with the lo- 
cal density approximation to extract the thermodynamic 
quantities in an inhomogeneous system. Our results will 
contribute to the worldwide experimental efforts in sim- 
ulating theoretical models of interacting systems using 
cold atoms in optical lattices. In this paper, we calcu- 
late local thermodynamic quantities such as local density, 
number of double occupation sites, entropy, and various 
compressibilities for the ID Fermi-Hubbard model. We 
present our results as a function of the chemical potential, 
particle number density and total number of atoms. The 
chemical potential and particle number density monoton- 
ically vary across the lattice due to the external harmonic 
potential so that the chemical potential can be considered 
as an effective spatial coordinate. At zero temperature 
and high temperatures, our work is complementary to 
some of the work presented in refs. [Ill4l7l |. However, 
we go beyond these studies by systematically extend- 
ing zero temperature calculations to any finite tempera- 
ture. We use both the Bethe ansatz and thermodynamic 
Bethe ansatz approaches. We also present an experi- 
mental scheme to deduce the calculated thermodynamic 
quantities from the column densities. Further, using the 
local density approximation, we calculate the entropy, 
compressibility, and total number of double occupation 
sites as a function of the total number of atoms in the 
lattice. 

The present paper is organized as follows. In section II, 
we present our model. In section HI, we present our zero 
temperature calculation procedure and the results. In 
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section IV, we discuss our finite temperature TBA calcu- 
lation and present various thermodynamic quantities for 
representative values of system parameters. In section V, 
we provide an experimental scheme to show the connec- 
tion between calculated thermodynamic quantities and 
experimentally measurable column densities followed by 
the calculation of thermodynamic quantities in a trapped 
system. Finally in section VI, we provide a discussion and 
draw our conclusions. 

II. THE MODEL 

We consider a two-component, cold Fermi gas in a ID 
optical lattice with N atoms and Ns lattice sites. We 
assume that the inter-site tunneling amplitude t is in- 
dependent of the pseudospin (cr). Here the pseudospin 
a =t,i refers to the two hyperfine states of the atoms. 
The system of fcrmions in the lattice can be described by 
the Fermi-Hubbard model: 



we consider an equal number of atoms in each pseudospin 
state. 




FIG. 1: Zero-temperature particle density as a function of the 
chemical potential /i at = 5. 



III. ZERO TEMPERATURE LIMIT 
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where c\^(cia) creates (destroys) a Fermi atom with 
pseudospin a at lattice site i. The density operator is 
^icr = c|^g.Cio-, < ij > indicates the nearest neighbor pair 
of sites and the chemical potential for fermions is /i. In 
general, the external trapping potential Vi = 7RI at site 
i for atoms can be treated in the local density approxi- 
mation by replacing /i by /i — . The on-site interaction 
U is sensitive to both laser intensity and the s-wave scat- 
tering length of the two Fermi species. As U is linearly 
proportional to the s-wave scattering length, it can be 
tuned to negative or positive values by exploiting Fesh- 
bach resonances. In the present study, we consider the 
repulsive Hubbard model so that U can have only positive 
values. In cold-gas experiments, the tunneling amplitude 
is controlled by the intensity of the standing laser waves. 
In current experimental setups, the tunneling amplitude 
is exponentially sensitive to the laser intensity while the 
on-site interaction is weakly sensitive. Thus, U jt can 
be controlled by a single parameter: the laser intensity. 
Further, standing laser waves allow for a variation of the 
dimensionality of the system. A reduced dimensional- 
ity is achieved by freezing the atomic motion in certain 
directions. Experimentally, the ID geometry can be re- 
alized by a strong confinement in the transverse direc- 
tion with an additional periodic potential applied along 
the other direction. The number of atoms in each pseu- 
dospin state is conserved because of the large energy gap 
between different spin states due to the external magnetic 
field present in most experiments. In the present study. 



The ID ground-state energy of the Fermi-Hubbard 
model can be calculated from a Bethe ansatz [s, 3, 19| . 
As a very good approximation, the ground-state energy 
per site E{n,U,t) — t x e{n,U/t — u) can be written 
based on the data obtained from a numerical solution of 
the Bethe- Ansatz [III, 111, H , 



2/(u) 



sm 
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This is valid in the region where < n < 1, where n = 
N /Ns is the average atom density per site. For the region 
1 < ri < 2, the energy can be obtained from the particle- 
hole transformation; e>(n,M) — [n— \)u + e<(2 — n,u). 
The function f{u) is determined from the equation 
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where Jm is the mth order Bessel function. We find the 
chemical potential by minimizing the grand canonical en- 
ergy under the fixed total particle numbers; fJ./t — de/dn, 



t 



2cos(^ 



< n < 1 



it-|-2cos( ^^^^ ), 1< n < 2. 



(4) 



The average particle densities are found by inverting this 
equation. 



3 



0, n/t < -2; 

^arccos[-M/2t], -2 < fi/t < -2 cos[7r//(u)] 

-2 cos[7r//(M)] <^i/t <u + 2 cos[7r//(u)] (5) 
2- ^arccos[(/^/i- u)/2], it + 2 cos[7r//(M)] < fi/t<u + 2 
2, n/t>u + 2. 
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FIG. 2: Zero-temperature compressibilities and double occu- 
pation as a function of chemical potential /i. The interaction 
is fixed at U /t = 5. 



For a representative value of the on-site interaction U /t = 
5, the zero temperature particle density as a function of 
the chemical potential is shown in FIG. [U As we stated 
before, the chemical potential-axis can be considered as a 
spatial coordinate where the minimum chemical potential 
is at the edge of the trap. As can be seen from the figure, 
five phases- two metallic phases, two insulating phases, 
and a vacuum phase- can coexist inside the trap if the to- 
tal number of particles in the trap is large. While density 
pins at 1 in the Mott insulating phase, it pins at and 2 
in the vacuum and band insulating phases, respectively. 

We calculate two compressibilities k* — dn/dfj, and 
K — K* /n. Panels (a) and (b) of FIG. [2]show these com- 
pressibilities as a function of chemical potential. The 
sharp peak of k at the edge of the trap is because 
of the low density. While this sharp peak is absent 
in K*, it shows two compressible regions correspond- 
ing to the lower metallic band and the upper metal- 
lic band. Further, we calculate the double occupation 
d = '^iini^Ui^) — de/dU and the double- occupancy com- 
pressibility Kc = dd/d^ which are shown in panels (c) 
and (d) of FIG. [H Notice that the panel (b) must be 
compared with the panel (d) as both represent the un- 
normalized compressibilities. Because of the sharp peak 
in the normalized compressibility k, we do not attempt 
to compare the normalized double- occupancy compress- 
ibility. Although we have fixed the on-site interaction 
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FIG. 3: Zero-temperature compressibilities as a function of 
the particle density for three representative values of interac- 
tions, long dashed (grey) line: U/t = 1, short dashed line: 
U/t = 5, solid line: U/t = 10. 



u to be 5, the qualitative behavior would be the same 
for all larger values of u. The three incompressible re- 
gions in panels (a), (b), and (d) corresponding to a pin- 
ning of the density at n =0, 1, and 2. Notice that the 
double- occupancy compressibility Kc measures only the 
core region [panel (d)] so that it is not sensitive to the 
edge of the cloud. Using a high temperature expansion 
Scarola et al. [31 have calculated the compressibility k 
and the double- occupancy compressibility Kc for a three 
dimensional optical lattice. In the next section, we show 
that our high temperature compressibilities in ID also 
show very similar characteristic behavior as those of the 
calculation in Ref . [ll] . In the insulating phases the com- 
pressibilities are zero so they are discontinuous at densi- 
ties n = 0, 1, and 2. These characteristic properties are 
directly related to the formation of gaps in the insulating 
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phases. As we will see later, these discontinuities smear 
out as the temperature increases. As shown in panel (d) 
of FIG.m the peak corresponding to the lower band dis- 
appears in Kc so that the double- occupancy compressibil- 
ity reveals the double occupancy and the Mott insulator 
transition in the center of the trap. The measurements of 
both double occupation and a quantity related to com- 
pressibilities, dd/dN — Kc/k* for a three dimensional 
lattice were made in Ref .ItI. In FIG. [3l we compare 
K* and Kc for different values of on-site interactions u. 
When the chemical potential is large, the upper band is 
occupied so that the system has doubly occupied sites. 

I 



For this identical to k* at chemical potentials 

corresponding to the upper band, but it is almost zero at 
the lower band. 



IV. FINITE TEMPERATURE RESULTS 

At finite temperatures, we numerically solve TBA 
equations to find compressibilities and entropy. Follow- 
ing Takahashi plj , the thermodynamic potential per site 
is given by 



n = Co — ^ — 



Po(fc)ln[I + e(fc)]dfc 
I 



ao(A)ln[l + 77i(A)]dA 



(6) 



The energy per site eo = 2i/, and two distribution func- 
tions of fc's and A's are given by [22[, 



po(fc) = hcos/c / ai(A - sinA:)CTo(A)(iA (7) 

1 

a-o(A) = — / s(A-sinfc)dfc. 



Here ai{x) — Au/['k{u'^ -f 16a;^)] and s{x) ~ 
csc(27rx/u)/u. The particle-hole ratios of k excitations 
and A excitations, ^(fc) and ?7i(fc) are obtained by an 
infinite set of nonlinear integral equations: 



In^(fc) 
In (A) 
In??; (A) 



Ko(fc) 

T 



j dAs{A - sink) ln(- 



! +??;( A) 



s*ln[l + 772(A)] - f s(A-sinA:)ln[l+^"l(/c)]cosfcd/c 

J — TT 

s*ln[l + 77^(A)] - / s(A-sin/s)ln[l+^(fc)]cosWfc 



(8) 
(9) 



and for j > 2, 



ln77,(A) = .*ln{[l + 77,_i(A)][l + 77,+i(A)]} 
ln77^.(A) = .*ln{[l + 77^_i(A)][l + 77;-+i(A)]}. 



(10) 
(11) 
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Here s*f{A) = /^^ s(A - A')/(A')dA' and Ko(fc) = 
~-2t cos fc - 4t dAs(A - sin k) x Rey^l - (A-ui/A)^. 
In order to calculate the thermodynamic potential nu- 
merically, one has to cut off the set of infinite equations 
at a finite number j . We achieve this by following the nu- 
merical procedure proposed by Takahashi et al. [23| . The 
infinite set of equations is truncated by replacing s(A) by 



5{A)/2 at j > ric- Then the integral equations are con- 
verted into a set of matrix equations in which 277^ + 1 
unknown functions are represented in terms of discrete 
points of k and A. These non linear matrix equations are 
then solved iteratively using Kepler's method. The de- 
tails of the numerical procedure can be found in Ref. [IS] . 
From the numerical solutions of the non-linear integral 



5 



2.0r 



s 1.0 



0.0^ 
-20 



A 

*j 

—' *n 


[C-^ 

M 


J' 

f =1 








2 2 4 6 



-10 



10 



20 



30 



FIG. 4: Particle Density as a function of the chemical poten- 
tial at U /t = 5. The four lines ranging from solid to short 
dashed represent different temperatures {kBT)/t = 0.1, 1, 2, 
and 10 respectively. 
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FIG. 5: The temperature dependence on entropy and par- 
ticle density at U/t = 5. Because of particle- hole symme- 
try, the entropy at n = 0.5(0.25) is equal to the entropy at 
n= 1.5(1.75). 



equations, we calculate thermodynamic potential fl, par- 
ticle density n = —dft/d^, compressibility k* = dn/d^, 
entropy S = —dil/dT, and double- occupancy compress- 
ibility Kc = dd/dfj,, where d = X]i("«t"4) ^ dfl/dU is 
the double occupation. 

First, we present particle density as a function of chem- 
ical potential at different temperatures. As shown in 
FIG. m the plateau region at half filling smears out as 
the temperature increases. However at very low temper- 
atures, the plateau region still survives showing a Mott 
insulating gap at n = 1. The calculated entropy at three 
different filling factors is shown in FIG.[S] In general, the 
entropy monotonically increases with temperature, how- 
ever, the increase is very strong at lower temperatures. In 
Ref . [i^l , the entropy at half filling for 3D optical lattice 
systems is calculated focusing on the anti-ferromagnetic 
(AFM) transition into Neel state. Even for a 2D finite 
lattice, a low temperature crossover scale is identified as 
a growth of AFM order 25|. As the finite temperature 
Neel transition is absent in ID systems and the system 
under consideration is in thermodynamic limit, we do not 
find a crossover scale in our ID calculation. 

The compressibility k* and the entropy S at fixed tem- 
peratures are shown in FIG.[S] As can be seen, the double 
peak structure in k* and S disappear at higher temper- 
atures. The disappearance of these sharp features may 
be used as a probe of temperature in experiments. The 
enhancements of compressibility contrast to the suppres- 
sions of entropy at lower temperatures are attributed to 
the metal-insulator transition at the center of the trap. 
At sufficiently high temperatures, the effect of interaction 
is weak. As the temperature decreases, the interaction 
leads to the opening of an Mott gap in the excitation 
spectrum at half filling. The low temperature results of 
FIG. m clearly indicate this characteristic feature. The 
wiggles in k* at ksT = O.lt are due to the numeri- 
cal noises in our low temperature calculation; however, 
the qualitative feature resembles the zero temperature re- 




FIG. 6: Compressibility and entropy at U/t — 5. The four 
lines ranging from dark to light colors represent different tem- 
peratures {kBT)/t = 0.1, 1, 2, and 10 respectively. 



suits. Notice that the entropy is largest in the vicinity of 
half filling at higher temperatures, it is larger in metallic 
phases at lower temperatures. This interesting feature 
opens the possibility of using entropy measurements to 
detect th metal-insulator transition at the center of the 
trap, as we see in the next section. 

In FIG. [71 we present the compressibility k* and en- 
tropy S* at a fixed temperature for different interactions 
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FIG. 7: Compressibility and entropy as a function of particle 
density at ksT = 0.5 t. The three lines ranging from dark to 
light colors represent different interactions U/t — 1, 5, and 10 
respectively. 



FIG. 8: Compressibility and double-occupancy compressibil- 
ity as a function of the chemical potential at fesT = 1.33 t. 
The three lines, solid, short dashed (black), and long dashed 
(gray) represent different interactions U/t = 10, 7.5, and 5 
respectively. 



involving different u. As expected, the compressibility 
decreases with increasing interaction and it is at a min- 
imum in insulating phases. The interaction significantly 
increases the entropy at the vicinity of the Mott insulat- 
ing phase (in the regions 0.5 < n < 1 and 1 < n < 1.5). 
However, deep in the Mott insulating phase, the inter- 
action has a very weak effect on the entropy. This is 
because the interaction increases the excitation gap and 
the density of low-energy states. For larger interactions, 
the excitation gap scales linearly with the interaction, 
but it is exponentially small for smaller interactions. 

The compressibility k* and double- occupancy com- 
pressibility shown in FIG. [51 As we have seen in 
zero temperature results, both k and almost iden- 
tical at chemical potentials corresponding to the upper 
metallic band, but Kc is zero at lower chemical potentials 
corresponding to the lower metallic band. A high tem- 
perature expansion method in three dimensions carried 
out by Scarola et al. [l^ also shows very similar quali- 
tative behavior as our ID exact results. As we discussed 
before, the occupation of the upper band is directly ob- 
servable through the measurements of double- occupancy 
compressibility. As can be seen from the figure, the upper 
band moves toward larger chemical potentials as the in- 
teraction increases. Again, this is because of the increase 



of the excitation gap with interaction. 



V. CONNECTIONS TO EXPERIMENTS 

This section is devoted to a discussion of the con- 
nection of our calculations to the experimental data. 
The density variations which can be measured using in 
situ images, can be used to deduce the bulk thermody- 
namic properties within the local density approximation 
(LDA). In order to study ID-lattice fermions in labora- 
tories, an array of ID-lattice geometries is formed with 
a 3D optical lattice. Apart from the underlying har- 
monic potential, the lattice potential has the form V = 
Vo±[cos'^{kx) -\- cos'^ (ky)] -\- Vbz cos^(fcz), with k — 2tt/X, 
where A is the optical-trap laser wavelength. For the 
system to be in a ID regime, the transverse laser inten- 
sity Vo± must be much larger than the axial laser inten- 
sity Vqz so that the transverse tunneling rate is negligi- 
ble compared to the axial tunneling rate. Further, the 
transverse confinement huj± >> kst, ep so that only the 
lowest transverse mode in each lattice tube is populated. 

The LDA assumes that the atomic system is locally 
homogeneous with local chemical potentials given by 
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Mi = Mo ~ ^i- Here /io is the chemical potential at 
the center of the trap. Then combining this chem- 
ical potential with the experimentally measured col- 
umn density profiles [n(y, z)]oi the ID-lattice bundles, 
the thermodynamic quantities can be extracted. The 
first step is the construction of the linear axial density 
n{z) = J n{y, z)dy. For non-lattice ID systems, these 
linear density profiles have been already measured at 
Rice University [2^. Notice that the compressibilities 
K and K* are related to the density n{z) and dn/d^. 
Since dn/d^ cx dn/dz, these compressibilities can be ex- 
tracted from the column densities. The extraction of the 
double- occupancy compressibility Kc = dd/dfj, involves 
the calculation of dd/dN. As presented in Ref. 0, this 
quantity can be extracted from the double occupancy. 
Then using the relations dd/dN — {dd/dfi){dfi/dN) 
and dN/dfi = {dN/dz){dfi/dz)^^, the double- occupancy 
compressibility can be extracted from experimental data. 

The entropy S = [dP/dT)^ needs two different con- 
figurations of pressure P at two different temperatures 
T, but at the same chemical potential m- Following 
the proposal by Ho and Zhou [27|, the local pressure 



P{z) = ma;^n(z)/(27r) can be extracted as for a 3D 
gas [2^. Two different column density profiles can be 
measured at two different axial trapping frequencies uj^i 
and uJz2, where ujz2 — i^zi- After taking an in situ im- 
age at ojzi, and then adiabatically changing ojzi — > i^z2, a 
second image can be taken. During this adiabatic change 
in the trapping frequency, both temperature and chem- 
ical potential will slightly change. Then calculating the 
pressure at two different positions such that /i^i = fiz2, 
the entropy S{z) = [P(z2) - P{zl)]/{T2 - Ti) can be 
extracted, where Ti and T2 are the temperatures of the 
system at the axial trapping frequencies uJzi and ujz2 re- 
spectively (27l.[28|. 

So far the present discussion has been about how to 
extract the homogeneous thermodynamic quantities from 
the experimental measurements for tapped atoms. Let us 
now turn to the calculation of the entropy and the com- 
pressibility in harmonic traps so that our calculation can 
be compared with the experimentally determined total 
entropy and compressibility of trapped Fermi gases. Con- 
sider the total number of atoms N — J dzn{z). Within 
the local density approximation, the integration over z 
can be converted to an integration over the local chemi- 
cal potential d/i = —mujzdz and we find 




FIG. 9: [Color onlinejTotal entropy as a function of the total 
number of atoms in a harmonic trap. 
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Figures [9] and [TO] respectively show the entropy S = 
St \/ muP' jit and k* = ymaP/^i as functions of the 
total number of particles N = N^Jmui'^ /2t for different 
interactions U /t &t a. fixed temperature fc^T = 0.5. Both 
the entropy and compressibility clearly show the appear- 
ance of insulating and metallic phases in the trap as we 
increase the number of atoms or the interaction strength. 
As an example, we plot the density profiles as a function 
oiz = muj'^z/2t in FIG.[n]for the case oiU ^ lOt. These 
density profiles correspond to different numbers of atoms 
N . Notice that each maximum of the compressibility in- 
dicates the existence of two metallic phases (upper band 
and lower band) while each minimum indicates the exis- 
tence of two insulating phases (Mott insulator and band 
insulator) inside the cloud. For example, when = 1, 
the compressibility does not show any minima so that the 
entire cloud is metallic. For the case of A^ = 8, the com- 
pressibility shows two maxima corresponding to the two 
metallic phases separated by the Mott insulator phase. 
Though it is very weak, the same behavior can be seen 
in the entropy profile too. 

The total double occupation Dt — J dzd{z) has the form 



M.J^I'-^r,. ,12) 

Here a; is the one-dimensional trapping frequency and 
ft = fi/t. The central chemical potential /io = Mo/^ is de- 
termined by the total number of atoms in the trap. Using 
the same change of variable, we find the total entropy and 
the total compressibility. 



2 J^ dfx. (14) 

The zero temperature total double occupation D = 
Dx^/muj'^ /2t as a function of N for different interac- 
tions is shown in FIG. [121 As one expects, the double 
occupation increases with the number of atoms but it is 
suppressed with increasing interaction. 
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FIG. 10: [Color onlincjTotal compressibility as a function of 
the total number of atoms in a harmonic trap. 
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FIG. 11: Density profiles for different atom numbers in the 
ID lattice for U/t = 10. 

VI. DISCUSSION AND CONCLUSIONS 

We have calculated several thermodynamic quantities 
for ID fermions in optical lattices. Modeling the fermions 
in the lattice by the ID Hubbard model, we have calcu- 
lated the particle density, compressibility, and entropy 
by means of the Bethe ansatz method. We find that 
these quantities have characteristic features that depend 
on the chemical potential, temperature, and interaction. 
As one expects, the sharp features on these quantities at 



low temperatures become smooth at higher temperatures 
and the visible signs of the Mott plateaus disappear. Our 
compressibility calculation shows that one can use com- 
pressibility measurements to detect stable compressible 
phases. In particular, we find that interaction effects, 
double occupation effects and temperature effects can be 
probed through measurements of the double- occupancy 
compressibility without the effect of trap edges. At rel- 
atively higher temperatures, while the compressibility 
K* — dn/d^i is finite in the entire range of the atomic 
cloud, the double- occupancy compressibility Kc = dd/d^ 
is finite only in the upper band. Therefore, the double- 
occupancy compressibility can be used to detect the oc- 
cupation of the upper band at the center of the trap. Fur- 
ther, we find that entropy also shows characteristic fea- 
tures; at higher temperatures entropy is larger at half fill- 
ing, at low temperatures it is larger at the lower and the 
upper bands. Therefore, entropy measurements would 
allow one to monitor the metal-insulator transition at 
the center of the trap. We use the local density approxi- 
mation to take into account the effect of the underlying 
harmonic trapping potential and find that measurements 
of entropy and compressibility of the entire cloud reveal 
the detailed structure (namely the existence of insulat- 
ing and/or metallic phases) of the trapped atoms in the 
lattice. 

In current cold-atom experiments, entropy can be 
monitored more easily than temperature. The lowest 
achieved entropy per particle is S ~ In 2. Our calcu- 
lations (FIG. [7]) show that entropy at half filing is ap- 
proximately equal to In 2. Therefore, the current lowest 
temperature experiments should allow one to observe the 
insulating gap and the metal-insulator transition. One 
dimensional lattice fermions with an attractive interac- 
tion are analogous to the half-filling repulsive Hubbard 
model with a finite magnetization. Therefore, the present 
work can be generalized to an attractive Hubbard system 
by introducing a population imbalance. 
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